Critical properties of the three-dimensional equivalent-neighbor model 
and crossover scaling in finite systems 
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, Accurate numerical results are presented for the three-dimensional equivalent-neighbor model on 

D ■ a cubic lattice, for twelve different interaction ranges (coordination number between 18 and 250). 

P_l ' These results allow the determination of the range dependences of the critical temperature and 

various critical amplitudes, which are compared to renormalization-group predictions. In addition, 
' the analysis yields an estimate for the interaction range at which the leading corrections to scaling 

vanish for the spin-i model and confirms earlier conclusions that the leading Wegner correction 
must be negative for the three-dimensional (nearest-neighbor) Ising model. By complementing 
' these results with Monte Carlo data for systems with coordination numbers as large as 52514, the 

full finite-size crossover curves between classical and Ising-like behavior are obtained as a function of 
a generalized Ginzburg parameter. Also the crossover function for the effective magnetic exponent 
is determined. 
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I. INTRODUCTION 



Over the past decades, several techniques have been applied to investigate how the critical behavior of systems 
depends on the range of the interactions. Before the general acceptance of the concept of universality, it was not 
at all clear that the critical properties of all systems with a one-component order parameter and ferromagnetic (i.e., 
attractive) interactions with a finite range are described by the Ising universality class. Since it was realized that most 
interactions in nature are not necessarily restricted to the nearest neighbors, one thus tried to determine the properties 
of models with a larger coordination number q. Another motivation, which plays a more important role in the present 
, work, is the fact that in the limit of infinite interaction range one recovers the classical or mean-field model. Since 
' the latter model can be solved analytically whereas no exact solution has been found for three-dimensional systems 
with a finite interaction range R, it is of interest to see how the crossover takes place from finite to infinite R. A 
natural choice for the examination of this crossover is the so-called "equivalent-neighbor" model, introduced by Domb 
and Dalton Q . In this generalization of the Ising model, each spin interacts equally strongly with all its neighbors 
within a certain distance, whereas all remaining interactions are equal to zero. In Ref. |l[, series expansions have 
been used to investigate two- and three-dimensional systems with interactions extending up to the third shell. On a 
simple cubic lattice this corresponds to 26 neighbors and on a face-centered cubic lattice even to 42 neighbors. While 
a general trend toward the mean-field properties, especially for the critical temperature, is clearly visible from these 
results, several problems emerge. First, with increasing interaction range, increasingly longer series are required to 
achieve a certain degree of convergence. Secondly, it appears that the maximum coordination numbers examined by 
this method are not large enough to observe the asymptotic deviations from the mean- field behavior Q. Although 
Ref. Q was published over 30 years ago, it appears that, especially in three dimensions, no substantial progress 
toward larger coordination numbers has been pursued. This is probably caused by the fact that also other techniques 
are plagued by serious difficulties upon increase of the interaction range. For example, Monte Carlo (MC) methods 
in general suffer from a serious decrease in efficiency if the number of interactions increases. Mon and Binder || 
studied two-dimensional (2D) spin systems with a maximum coordination number q = 80, compared to q = 12 and 
q = 18 for quadratic and triangular lattices, respectively, in Ref. |Q. Furthermore, they derived the R dependence 
of critical amplitudes from scaling considerations. However, it still proved difficult to reach the asymptotic regime 
where the predictions were expected to hold. In a subsequent paper H], Luijten et al. confirmed the predictions of 
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Ref. [|| from a renormalization-group (RG) analysis and revealed the existence of a logarithmic R dependence in the 
shift of the critical temperature. Thanks to the advent of a dedicated MC algorithm for long-range interactions [0, 
systems with large coordination numbers could be simulated without loss of efficiency. Thus, in the same paper the 
critical properties were determined for quadratic systems with coordination numbers up to q = 436. It was explicitly 
verified that all examined systems belong to the 2D Ising universality class and the predicted R dependence of the 
critical amplitudes could indeed be observed, as well as the approach of the critical temperature toward its mean-field 
value. It is the purpose of the present work to extend this analysis to three-dimensional (3D) systems. Apart from the 
possibility to verify the predicted range dependences in three dimensions, a precise knowledge of the critical properties 
of spin models with an extended range of interaction also serves a further purpose. Namely, it allows the study of two 
forms of crossover in these systems. Finite-size crossover only pertains to finite systems at the critical temperature 
and denotes the transition from the classical regime where the interaction range is at least of the order of (some 
power of) the system size to the nonclassical (Ising) regime where the system size is much larger than the interaction 
range. Thermal crossover, on the other hand, occurs when the temperature is moved away from its critical value. 
The interplay between the range R of the interactions and the decreasing correlation length £ determines the location 
of the crossover to classical critical behavior. If R is small, the temperature distance to the critical temperature T c 
must be made rather large before £ and R are of the same order of magnitude. In such systems, no crossover to 
mean-field-like critical behavior can be seen, because one has already left the critical region. However, for R large, it 
is well possible to observe both Ising-like and classical critical behavior. This dependence on both t = (T —T c )/T c and 
R is expressed by the Ginzburg criterion . Both variants of crossover have been studied for 2D systems in Refs. > 
which showed that accurate information on crossover scaling functions can be obtained by numerical techniques. In 
the light of a comparison to experimental results on the one hand and theoretical calculations of crossover scaling 
functions on the other hand, it is extremely relevant to investigate the 3D case as well. Here, I present the results of 
MC simulations of systems with interactions up to a distance of \/l4 lattice units (thirteenth shell), which corresponds 
to 250 equivalent neighbors. Although larger interaction ranges do not diminish the efficiency of the MC algorithm, 
an accurate determination of the critical properties for larger R is hampered by a different effect. Indeed, such a 
determination is only possible in the Ising limit, which implies that the smallest linear system sizes must be of the 
order of L min = 0(R i ^ 4 ^ d ' ) ) 0], where d indicates the dimensionality. Thus, for d = 3 the smallest allowable systems 
contain of the order of R 12 spins and one can only hope that this relation exhibits a prefactor considerably smaller 
than unity. 

The results of the MC simulations arc then used to determine the finite-size crossover functions for several quantities. 
It should be noted that for a full mapping of these function very large coordination numbers are required: simulations 
have been carried out for q up to 52514. Yet, an independent determination of the critical temperature of these 
systems is not required, but can be obtained by extrapolation. It suffices thus to study modest (L < 40) system 
sizes for these interaction ranges. The determination of thermal crossover functions will be the subject of a future 
paper || , as it requires calculations which are actually complementary to those of the present work (results for the 
susceptibility can be found in Ref. jLO)). Indeed, for a determination of the critical properties by finite-size scaling 
and for the mapping of the finite-size crossover functions, all data must lie within the finite-size regime, whereas for 
thermal crossover scaling care must be taken that the data lie outside this regime. 

Two further questions that are addressed in this paper concern the corrections to scaling. In the first place, the 
range dependence of the thermal finite-size corrections is shown to be in very good agreement with the predictions of 
Ref. |j . Secondly, the finite -size corrections due to the leading irrelevant field are analyzed and the related variation 
of the (j) 4 coefficient in the Landau-Ginzburg- Wilson (LGW) Hamiltonian is obtained. This permits an estimation of 
the interaction range for which this coefficient coincides with its fixed-point value and confirms that for the three- 
dimensional nearest-neighbor Ising model it does not lie between the Gaussian fixed point and the Ising fixed point. 

The outline of this paper is as follows. In Sec. O, I briefly summarize earlier predictions for the range dependence 



of critical amplitudes and discuss the shift of the critical temperature as a function of interaction range. Section III 
gives details of the Monte Carlo simulations. Furthermore, the determination of the critical temperatures is discussed 
as well as the analysis of the range dependence of corrections to scaling. The variation of critical amplitudes as a 
function of interaction range is treated in Sec. ^ and finite-size crossover curves are obtained in Sec. |v|. I end with 



some concluding remarks in Sec. VI 



II. SUMMARY OF RENORMALIZATION-GROUP PREDICTIONS 

In the absence of an external field, the equivalent-neighbor or medium-range model is defined by the following 
Hamiltonian, 
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H/k B T = -Y,K{v l - r 3 )s lSj , (1) 

m 

where s = ±1, the sum runs over all spin pairs and the spin-spin coupling is defined as K(r) — J > for |r| < R m 
and K(r) = for |r| > R m - I first summarize the findings of Ref. Q for the R dependence of critical properties, as 
obtained by an RG analysis. Although at first sight this approach is not very different from a simple scaling analysis, 
it offers several advantages. The formulation in terms of two competing fixed points provides a clear insight into the 
crossover mechanism: for R large the coefficient of the tfi 4 term in the LGW Hamiltonian is suppressed with respect to 
the quadratic term in this expression. Thus, the renormalization trajectory passes close to the Gaussian fixed point 
and the critical amplitudes pick up a specific R dependence which is determined by the flow near this fixed point. For 
any finite R, the system will still flow to the neighborhood of the nontrivial (Ising) fixed point (cf. Fig. 1 in Ref. [Q). 
However, the R dependence reveals some aspects of the Gaussian fixed point which are not normally seen in Ising-like 
systems. For example, near this fixed point the thermal exponent yt and the leading irrelevant exponent yi assume 
the values 2 and 4 — d, respectively, which coincide for d — 2. Such a coincidence would lead to logarithmic factors in 
the scaling functions, were it not that the Gaussian fixed point is unstable for d — 2. In contrast, the R dependence of 
scaling functions indeed allows the observation of such logarithms. The occurrence of these dependences is not easily 
found from a scaling analysis. 

For the magnetization density m and the magnetic susceptibility x the following range dependences have been 
obtained, 

m oc t P R (^-d)/{A-d) f (2) 

X OC t -7flM(l-7)/(4-«0 ; (3) 

where (3 and 7 denote the standard Ising critical exponents. Furthermore, the finite-size scaling behavior of these 
quantities was derived as 

m = L y h -d R (M-4y h )/(4-d) £(1) ^ R -2(2y t -d)/(4-d) R -*m/(4-*> jhL VK R (3d-iy h )/(4-d)^ ^ (4) 

X = tfyh-d R 2CAd~4y h )/(4-d) j(2) / f £j/ t R -2(2y t -d)/(4-d) ^ RVi R -4y, /(4-rf) hL y h R (3d~4y h ) / (4-d) \ _ ^ 

Here, /J^ denote universal scaling functions, yt and are the thermal and leading irrelevant exponent introduced 
before, and y^ is the magnetic exponent, u and h are the irrelevant and the magnetic scaling field, respectively. 

Also the shift of the critical temperature with respect to its mean- field value has been calculated in Ref. ||. However, 
this treatment left several questions unanswered, which I will consider here in some more detail. A clear understanding 
of the nature of this shift is of particular significance for the crossover scaling, since one has to calculate the critical 
temperatures for systems with large coordination numbers by means of extrapolation. It was derived that under a 
renormalization transformation the contribution of the <fi 4 term to the quadratic term in the LGW Hamiltonian leads 
to a range-dependent shift of the reduced temperature t = (T — T c )/T c . For d = 2 it was found in Ref. || that this 
shift has the form 

T c -Tr= C ° + * lnR + - ■ (6) 
R z 

where Co and ci are constants. This expression has also been confirmed numerically, see Fig. 4 in Ref. Q. Interestingly, 
this result has been recovered in Ref. plf , where in addition it was found that the constant c\ has a universal 
value —2/ir w —0.6366. Indeed, this agrees with the value —0.624 (7) obtained from an analysis of the available 
data for 25 < R 2 < 70. (The somewhat lower value 0.609, corresponding to the coefficients quoted in Ref. ||, can 
be explained from the influence of the data point at R 2 — 16.2.) However, the result for general 2 < d < 4, a shift 
proportional R- 2d /( 4 ~ d \ clearly contradicts the results obtained from systematic expansions in terms of the inverse 
coordination number (but see the remarks at the end of this section). Brout p2| obtained to leading order a shift 
of the form 1/q oc 1/R d . This result was recovered by Vaks et al. |l3] and Dalton and Domb [Q. As indicated in 
Ref. [fijl, such an additional and actually dominant shift can also be obtained from the RG analysis by allowing for 
a (spherically symmetric) lower-distance cutoff a in the spin-spin coupling -K'(r). In momentum space the coupling 
then takes the form 

*(k) = c(-) J d/2 (kR)-c(-) (-) J d/2 (ka), (7) 

where c = JR d and J v is a Bessel function of the first kind of order v (cf. Eq. (A3) in Ref. Q|). The second 
term in this expression yields an additional contribution to the quadratic term in the LGW Hamiltonian, which is 



3 



precisely the 1/R d shift obtained by Brout. Furthermore, it contributes to the fc-dependent part of this term, which 
via the rescaling of the field (see Ref. |Q) leads to a 1/R d+2 shift. Note that, upon expansion in powers of R, a 
formulation in terms of the coordination number q implies such a shift as well. At even higher order, one finds (at 
rational dimensionalities) additional lni? dependences, as was first recog nized by Thouless [§.[] He has studied a 
modified form of the Ising model, where the system is divided into cells within which the spin-spin interactions are 
constant. The shift of the critical temperature as a function of the cell size is then calculated by means of perturbation 
theory. In three dimensions, the leading-order result of Brout is recovered, namely a shift proportional to l/q. In the 
next-to-leading term a logarithmic dependence on the coordination number is obtained, 

q q- 

whereas for d = 2 the logarithm emerges already in the leading term, 

The latter expression is in perfect agreement with Eq. (||) , whereas the logarithm in the higher-order term in (^J) has 
not been found in Refs. |^,^5|. Since the logarithms in Eqs. (||) and (^|) apparently follow from the same mechanism and 
the factor In R in Eq. (^J) is specific for the two-dimensional case (where all higher-order terms in the LGW Hamiltonian 
are equally relevant), I conclude that there must be two different sources for the logarithms, which happen to yield 
the same effect in d = 2. Indeed, the logarithms in Eqs. (||) and (||) arise from counter terms canceling the infrared 
divergences in the perturbation expansion. This appears to be intimately linked to the infrared divergences occurring 
in massless super-renormalizable field theories at rational dimensionalities Actually, the treatment of Ref. Q 

does account for logarithmic factors in d = 3, although at much higher order. For systems with a large interaction 
range, the first part of the renormalization trajectory passes close to the Gaussian fixed point. Near this fixed point, 
only the (j> A term is relevant for d = 3 and all terms (f> n with n > 6 are irrelevant. The marginal character of the 4> 6 
term produces a logarithmic range dependence in the shift of the critical temperature. However, since this logarithm 
stems from the term quadratic in <fi 6 and the field <fi is rescaled by a factor R , this contribution is extremely weak. 
An actual calculation shows that it leads to a shift proportional to In R/R 18 (x\nq/q 6 . In addition, the 4> 6 term will 
yield a correction of order R~ 8 . However, it may be added that it is generally expected that such high composite 
operators have very little influence near the Ising fixed point. 

Let me now briefly return to the leading shift R~ 2d /^ d ) obtained in Ref. It is instructive to note that this shift 
is consistent with crossover arguments first given by Riedel and Wegner jl8|]. Indeed, the Ginzburg criterion states that 
a crossover from classical to nonclassical critical behavior occurs as a function of the crossover parameter t^ 4 ~ d ^ 2 R d . 
In terms of a more general formulation, this parameter is written as t^/g, with <fi — (4 — d)/2 and g = R~ d . The 
crossover exponent 4> (not to be confused with the field <fi), introduced in Ref. Jig] , is just the exponent 4 — d 
of the relevant operator driving the system away from the Gaussian fixed point (i.e., the </> 4 term in the LGW 
Hamiltonian), divided by the thermal exponent yt = 2. Then, on general grounds pJL[L9|, the shift of T c is predicted 
to scale as g 1 ^ oc R 2d /( 4 d \ This is another indication that the shift terms in Eq. (H) originate from a different, 
complementary mechanism. In addition, it is noted that the formulation in terms of the crossover exponent <fi can be 
carried even further (see, e.g., Ref. ]20|]). Indeed, for any thermodynamic quantity P which is near the Ising critical 
point proportional to t Xl , the combined dependence on g and t will be 

Paj^- 1 ')^ 11 , (10) 

where xq denotes the t dependence of P near the Gaussian fixed point. In terms of t and R, this can be written as 

p K R d( Xl -xa)/^ t xi 5 (11) 

which yields, e.g., m cx #2d(/3-i/2)/(4-d) t /3 and x a i? 2rf(i- 7 )/(4-rf) t - 7j recovering Eqs. (§) and (|). 

III. MONTE CARLO SIMULATIONS 



1 This work only came to the attention of the author after the publication of Ref. 
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A. Simulational details 



I have carried out extensive simulations of 3D simple cubic lattices consisting of L x L x L lattice sites with periodic 
boundary conditions. Each spins interacts equally with its q neighbors lying within a distance R m , i.e., the system is 
described by the Hamiltonian (|l|). For the simulations I have used the cluster MC algorithm introduced in Ref. ||. 
Its application to the present case is described in more detail in the appendix of Ref. 0. In order to avoid lattice 
effects I formulate the analysis in terms of an effective interaction range R j| , 

r2 _ E^(r*-r 3 ) 2 ^- 

= -$^|ri-r 3 -| a with \n - rj\ < R m . (12) 

It is easily seen that lim^^oo R 2 — 3i?^/5. Table | lists R m , q, and R for the first thirteen neighbor shells which have 
been examined in the present work. 

Several tests have been carried out to check the implementation of the algorithm. For R 2 n — 1 exact results (for 
L = 3,4) and accurate MC data are given in Ref. |2l| and for R 2 n = 2,3 alternative MC programs were available, 
allowing the verification of the data for various system sizes. I have carried out very long Monte Carlo simulations 
(10 9 and 10 s Wolff clusters, respectively) for L = 4 and L — 20 for these ranges, at couplings close to K C (R). On 
the other hand, if one takes into account all lattice symmetries, an explicit summation over all states is feasible for 
L = 3 (2 27 sa 1.34 x 10 8 configurations). For this case, I have carried out simulations for all ranges 1 < R 2 m < 14. No 
systematic deviations could be observed. The actual simulations were carried out for systems up to L = 200 (8 million 
spins); the number of samples was chosen depending on the system size. As a rule of thumb, the amplitude ratio Q 
(to be defined below) had a relative accuracy on the permille level for the largest systems. 

B. Determination of the critical temperatures 

In order to analyze the range dependence of several quantities, an accurate knowledge of the critical temperature for 
each single value of R m is required. The critical temperatures of systems with interaction ranges corresponding to the 
first thirteen neighbor shells have been determined using the amplitude ratio Ql — ( to2 )|/(w 4 )l. For the 3D Ising 
universality class and a cubic geometry with periodic boundary conditions, this quantity has, in the thermodynamic 
limit, the universal critical-point value Q = Qi = 0.6233 (4) pl[ . As mentioned in Sec. |, an accurate determination 
of the critical point is mainly hampered by the requirement that one must reach the Ising limit, i.e., L min ss i? 4 . 
For the inner shells, the smallest system sizes that could be used in the finite-size analysis were of the same order as 
in an analysis of the 3D nearest-neighbor Ising model, i.e., L > 5. For the remaining shells, the smallest allowable 
system sizes, as determined from the quality of the least-squares fits, followed the restriction L > i? 4 rather closely. 
Only for the outermost shells this criterion could be slightly relaxed. Thus, the accuracy of the fit results decreases 
considerably with increasing interaction range, because the finite-size data cover a much smaller range of system sizes 
and all the accurate results for small system sizes have to be excluded from the analysis. The least-squares fits were 
made using the finite-size expansion for Q given in Ref. |2l[| , 

Ql(K) = Q + ai(K - K C )V* + a 2 {K - K c ) 2 L 2yt + ■■■ + biU" + b 2 L 2 ^ + ■ • • , (13) 

where K denotes the spin-spin coupling, K c the critical coupling, and the dj and hi are nonuniversal (range-dependent) 
coefficients. The exponents yt and yi are the thermal and leading irrelevant exponent, respectively. They are approx- 
imately given by: yt = 1.587 (2) and yi — —0.82 (6) |piJl , where the latter exponent was kept fixed in all analyses. 
Table ffl shows my resulting estimates for Q and K c . In the first place, one notes that all systems belong, within the 
statistical accuracy, to the 3D Ising universality class. The critical couplings for the first three shells are in agreement 
with the old series-expansion results of Domb and Dalton. In order to improve the accuracy of the results, I have 
repeated all analyses with Q fixed. 

The results of the finite-size analyses permit some additional tests of the scaling predictions of Refs. Indeed, 
the range dependence of the thermal coefficient a\ in Eq. ([l3]) should take the same form as the first argument of the 
universal scaling functions (|j) and (0) . Upon expansion of such a scaling function one finds a temperature-dependent 
argument of the form atL^R- 2 ^^^ 4 -^ » -a[(K - K c ) / 'K c ]L yt R-^vt-d)/(A-d) ^ where a is a constant tna t does 
not depend on R. Thus a x = -aR- 2 ^i.-d)l^-i) jK c _ R -2(2 Vt -d)/(4-d) x R d _ ^2.652 _ g o.884_ Figure [I] shows a x as 
a function of the coordination number q. Both a curve ~ ^0.884 anc j a re f erence j me w itfi slope 1 are shown; evidently 
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the former describes the numerical data very well. Deviations for relatively small q are not disturbing, since the RG 
predictions are only valid in the limit of large interaction ranges and the small-q data may also exhibit some lattice 
effects. 

Of particular interest is also the range-dependence of the coefficient b\ in Eq. ([l3|), because this coefficient is 
proportional to (u — u*)/u*, where u is the coefficient of the 4> A term in the LGW Hamiltonian and u* is its fixed- 
point value Iffi] . As such, b\ yields information on the R dependence of the size and sign of the corrections to 
scaling that appear in the Wegner expansion p3| . This expansion describes the singular corrections to the asymptotic 
temperature dependence of thermodynamic quantities close to the critical point. For example, if u/u* > 1 the leading 
coefficient in the expansion for the susceptibility will have a negative sign and hence the susceptibility exponent 7 will 
approach the fsing value from above, cf . also Ref . [^4j . On the other hand, if u lies between the Gaussian and the Ising 
fixed point, i.e., < u/u* < f, the sign of the first Wegner correction will be positive and 7 will approach the Ising 
value from below. In order to extract the R dependence of u from the coefficient b\, the RG scenario of Ref. Q has to 
be reconsidered. It can be shown that in the large- R limit u = uq/R" 1 . Because uq will exhibit a remaining, weak R 
dependence for small R, I write it as Uo(R). The first part of the RG transformation is just a scale transformation in the 
neighborhood of the Gaussian fixed point, which cancels the factor 1/ R A in u. The <j) A coefficient can now be written as 
uq(R) = u* + [uo(R) — u*], which close to the Ising fixed point scales as u — > u' = u* + [u (R) - u*]L Vl R~ 4yi/ '( 4 ~ d ) [Q. 
Thus, the coefficient 61 in Eq. (|) is equal to c[u(R) - l]_R- 4 w/( 4 - d ), where u(R) = u (R)/ u* and c is a nonuniversal 
proportionality constant. For R large, u(R) should go to a finite constant and hence b\ is expected to be proportional 
to R- 4 Vi/( 4 ~ d ) in this limit. Just as for most other quantities, it is difficult to accurately determine b\ for large 
interaction ranges, because the small system sizes have to be omitted from the analysis. Nevertheless, the results 
shown in Fig. g appear to be well compatible with the predicted R dependence, with c[w(oo) — 1] w —0.14 (the latter 
estimate relies on the assumption that the asymptotic limit has actually been reached for the largest ranges shown in 
the figure). Unfortunately, no estimate for uq{R) for either R = 1 (nearest-neighbor Ising model) or any other R is 
known to the author, so that the overall constant c [which would have permitted the calculation of uq(R) from b± (R)] 
cannot be determined (cf. also Ref. p5|]). 

On the other hand, an estimate of the interaction range where u(R) = 1 does not depend on c, and so it can 
be predicted with a reasonable accuracy that this condition is fulfilled at R 2 ss 1.56. The interest of this point 
lies in the fact that the leading corrections to scaling should vanish there, which in principle allows a much more 
accurate determination of critical properties from numerical simulations. This approach was used for the first time in 
Ref. pH , where, amongst others, a spin-i model with nearest-neighbor coupling K nn and third-neighbor coupling K 3n 
was simulated. The ratio K% n /K nn was set to 0.4, which in hindsight proved to be somewhat too strong for fully 
suppressing the leading corrections to scaling. A newer estimate yielded K^/K nYl = 0.25 (2) as an optimal choice p6[ . 
Further studies of these systems were presented in Ref. [M , where the coupling constant ratio was systematically 
varied in order to eliminate the leading finite-size corrections. This lead to an estimate of K^ n /K nTl « 0.27. Both 
estimates turn out to be in quite good agreement with my prediction for general interaction profiles. Indeed, as follows 
from Eq. (12), an effective interaction range R 2 = 1 .56 can be obtained by, e.g., nearest-neighbor and next-nearest 
neighbor interactions with K2 n / K nn — 0.64 or by nearest and third-neighbor interactions with K3 n /K nn = 0.29. This 
also explains the finding of Ref. |^f| that Kin had to be chosen much larger than K^ n to reach the same effect. 

In this context it is of interest to review some series-expansion results for the leading correction amplitudes for 
the magnetization, the susceptibility, and the correlation length on simple cubic (sc), body-centered cubic (bec), and 
face-centered cubic (fee) lattices. Liu and Fisher [ p2[ concluded that the leading correction amplitudes are negative 
for the sc and bec lattices and gave various arguments that this also holds for the fee lattice. Furthermore, they argue 
that these amplitudes should vanish monotonically with coordination number (q — 6,8,12, respectively). This is 
indeed confirmed by the fact that the data in Fig. monotonically approach the predicted asymptotic R dependence, 
apart from statistical scatter. However, from the fact that for the sc lattice with q = 18 (R 2 — 5/3) the finite-size 
corrections have already changed sign, it would be expected that the correction amplitudes for the fee lattice are 
close to zero. In contrast, both the results of George and Rehr |^] and Liu and Fisher [^9), see Table III, exhibit 
a relatively weak variation with coordination number. On the basis of these results one would certainly expect the 
leading corrections to vanish at much higher coordination numbers. Thus, I conclude that, apart from the dependence 
on q (or R), the value of u has a rather strong dependence on the lattice structure as well. For completeness, it may 
be remarked that the analyses of the Monte Carlo data for the magnetization density and the susceptibility have 
revealed the same monotonic R dependence of the leading correction amplitude as that of the quantity b\ discussed 
above. 



IV. RANGE DEPENDENCE AT CRITIC ALITY 
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A. Critical temperature 



The estimates for the critical coupling as given in Table |l| can in principle be used to verify the predictions for the 
shift of the critical temperature. Because lattice effects are still relatively strong for the interaction ranges studied 
here, the coordination number q, appearing in, e.g., Eq. (jS), cannot be used directly. It is expected that these lattice 
effects disappear when the effective interaction range R is used instead. Thus, the predicted shift is rewritten as: 

rn-l- v 1 , c ° i Cl i c 2 + c 3 InR 

T c =qKc = l + - + -+ R6 +■■■, (14) 

where I have used the inverse critical temperature to conform to the earlier literature. Unfortunately, it turns out that 
even in terms of R the numerical data exhibit remarkably strong scatter for R^ < 10, making it impossible to obtain 
a sensible fit for the smaller interaction ranges. On the other hand, for R^ > 10, Eq. (|l4| ) describes the data very 
well. Because of the small variation of the In R term over the fit range, it was not possible to discern the coefficients 
C2 and C3. Thus, I have omitted c 2 altogether, which implies that this coefficient is absorbed into an effective value 
of C3. The resulting fit yielded the values Co = 0.498 (2), c\ = —5.7 (7), and C3 = 7.1 ( 9). Clearly, the last two 
estimates suffer from the fact that (for the available values of R) the last two terms in Eq. (jlj) lie quite close. Thus, 
it cannot be excluded that the high values of c\ and C3 are partially caused by a mutual cancellation and that apart 
from the quoted statistical errors there is a considerable systematic error. Nevertheless, as will be seen below, the 
accuracy of the resulting expression is sufficient to obtain rather precise estimates for systems with larger interaction 
ranges. In fact, if the results for R^ — 9, 10 are also included in the least-squares fit and the lattice effects are simply 
ignored, an essentially phenomenological interpolation formula is obtained, which for larger ranges turns out to agree 
very well with the first fit. 

In Refs. H,B|, series-expansion estimates are given for the coefficients Co and C3 in (14). In terms of an expansion 
in q, Dalton and Domb found the value 4.46 for the leading coefficient (confusingly, in later work jl],[30| the value 3.5 
was quoted) and for the prefactor of the logarithm Thouless obtained —2000/27 « —74.1. To compare these values 
to c and c 3 , I write ^R 3 m m §7r(f) 3 / 2 i? 3 » 9.013i? 3 . This yields c = 0.495 and c 3 = -2.74. In view of 

the various approximations that have been made, the agreement for c is truly remarkable. Because of the above- 
mentioned cancellation effects and because of the omission of c 2 in the fit, a sensible comparison for C3 is not possible. 
However, we note that also Thouless finds a relatively high value for c 3 . Figure || shows the various predictions for 
the shift of the inverse critical temperature. 



B. Magnetization density 

In the Monte Carlo simulations, I have sampled the absolute magnetization density (|m|). The dependence of this 
quantity on both L and R is given by Eq. (^) , from which the following finite-size expansion can be derived, 

m L (K, R) = V h ~ d {d Q (R) + di(R)[K - K c (R)}L Vt + d 2 (R)[K - K c (R)] 2 L 2yt + ■■■ + ei (R)L v ^ + • ■ •} . (15) 

For each single value of R, I have fitted the numerical data to this expression. The critical couplings obtained from 
this analysis are in agreement with those shown in Tabic ||. The corresponding estimates for yh are listed in Table IV 



The slight tendency of the estimates to decrease with increasing R, as well as the increasing uncertainties, can be 
explained from the requirement that the smallest system size included in the analysis must increase with R. When the 
analyses were repeated with the critical couplings fixed at the best values in Table [n| the agreement of the estimates 
for yh (also shown in Table IV) with the 3D Ising value yh = 2.4815 (15) was even better. Thus, this confirms 



the expectation that all these systems belong to the Ising universality class. The critical amplitudes do(R) can be 
used to extract the leading range dependence of the magnetization density. In order to maximize the accuracy in 



these amplitudes, the results shown in Tabic IV were obtained with the exponents yh and y t fixed at their Ising values 
(but the critical coupling K c was included as a free parameter). A fit of da(R) to the form do(R) = dR x for the 
largest three values of R yielded x — —0.87 (5), somewhat (although not significantly) smaller than the predicted 
value x = (3d — 4y/ l )/(4 — d) = —0.926 (6). This shows that the asymptotic regime, where higher-order corrections 
can be neglected, has not yet been reached. In general, the corrections are powers of R~ 2 [Q: 

do (R)=dR*(l+^ + ^ [ + --) (16) 



R 2 R 4 

Expression ( |l6| ) with one correction term allowed me to obtain a very acceptable fit (% 2 /DOF w 0.6) for all data 
points with 2 < R^ < 14 and yielded x = 0.923 (5), in excellent agreement with the RG prediction of Ref. Q. Figure || 
shows the MC results for do (R) together with the asymptotic range dependence and the full fit to the renormalization 
expression. 
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C. Susceptibility 



At criticality, the magnetic susceptibility is directly proportional to the average square magnetization. Thus, I have 
fitted the Monte Carlo data, for each interaction range separately, to an expression of the form 

X l(K, R) = s + L 2yh ~ d { Po (R) + Pl (R)[K - K C (R)]L* + p 2 (R)[K - K C {R)] 2 L^ + • • ■ + q 1 (R)L y ' +■••}, (17) 

where the additive constant sq originates from the analytic part of the free energy. In the further analysis, this 
constant has been set to zero, because it tends to interfere with the leading irrelevant term q\{R)L Vi . Just as for the 
absolute magnetization density, I list estimates for the magnetic exponent yn (Table |V|). Although, as expected, the 
uncertainty increases with R, one observes that all estimates agree with the Ising value. Also the critical couplings 
agree with those obtained from the fourth-order amplitude ratio and (|m|). Thus, I have repeated all analyses with 
K c fixed; the corresponding results for y^ are shown in Table ^ as well. Finally, I have fixed the magnetic exponent at 
yh = 2.4815 (but included K c as a free parameter) in order to obtain accurate estimates for po(R) (Table ^). Fitting 
a straight line pR~ x to the last three values yielded a slope —1.73 (9), which is consistent with the prediction —1.852 
[Eq. (0)]. A fit formula with one additional correction term, pR~ x (\ + bR~ 2 ), allowed the inclusion of several more 
data points and yielded x = —1.92 (11). Both fits and the numerical data are shown in Fig. ||. 

D. Connected susceptibility 

In principle, the connected susceptibility, given by 

k B T 

can be treated in the same way as the absolute magnetization density and the susceptibility. The main drawback 
of this quantity, being the difference of two fluctuating quantities, is that its statistical accuracy is relatively poor. 
Nevertheless, the magnetic exponents extracted from the numerical data for the individual interaction ranges are 
consistent with the Ising value and the finite-size amplitudes can be used to determine the range dependence of the 
connected susceptibility. As shown in Ref . Q , knowledge of this dependence is very useful to determine the thermal 
crossover curve for the susceptibility (which for T < T c is represented by x) from data for different R, because it 
makes it possible to divide out the subleading range dependence of this curve. Rather than giving the full details 
of the analysis, I restrict myself to Fig. ^, which shows the critical amplitudes together with the RG prediction 
fitted to it. Instead of x, the so-called scaled susceptibility ksTx is often considered. It has been noted for the 
two-dimensional case [9 , that the latter quantity exhibits considerably stronger deviations from the asymptotic range 
dependence, which are caused by the shift of the critical temperature. Figure ^| confirms that this also holds for the 
three-dimensional case. 

V. FINITE-SIZE CROSSOVER 
A. General considerations 

As stated in the introduction, the critical properties of the equivalent-neighbor models obtained in the previous 
section can now be used to find the finite-size crossover scaling functions describing the crossover from a finite mean- 
field- like system to a finite Ising-like system at T = T c . A detailed description of this phenomenon has been given 
in Ref. M. Qualitatively this crossover can simply be understood from the observation that systems with a linear 
size of the order of the interaction range are essentially mean-ficld-like systems, which are turned into systems with 
a short-range interaction if the system size grows beyond the interaction range. RG considerations have shown 
that the crossover is ruled by a generalized Ginzburg parameter G = LR^^'^~ d \ so that the mean-field regime 
corresponds to G <C 1 and the Ising regime to G > 1. The expression for G has also been obtained in Ref. [j3l| . 
It is numerically not feasible to observe the entire crossover regime in a system with fixed R by merely varying the 
system size, since it spans several decades in the parameter G. Thus, I construct the crossover curve by combining 
the results for various interaction ranges, just as this has been done for the two-dimensional case in Ref. || and for 
the three-dimensional thermal crossover in Ref. ]l0| . Since it turns out that for L < 20 the curves are affected by 
nonlinear finite-size effects, the smallest value of the crossover parameter that can be reached with the interaction 
ranges studied in the previous two sections is 20/(9. 168) 2 » 0.24. The true mean-field regime, however, is only reached 
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for much smaller G — 0(1CP 4 ). Thus, I have carried out simulations for systems with effective interaction ranges up 
to R 2 = 323.81 (R 2 n = 540), corresponding to coordination numbers as large as q — 52514. Evidently, the Monte 
Carlo algorithm introduced in Ref. || comes to its full glory here: The simulation of three-dimensional systems with 
so many interactions present would not have been feasible with either a Metropolis-type algorithm or a conventional 
cluster-building algorithm. The actual crossover curves shown below are obtained from a combination of the data 
for 2 < < 14, with system sizes between L = 20 and L = 200, and additional data for 20 different interaction 
ranges 18 < R^ < 540 . For the latter systems, the critical coupling was determined using the extrapolation formula 
discussed in Sec. IV A and subsequently simulations were carried out for 20 < L < 40 at each single value of Rf n . 

An additional complication is formed by the regime G> 1. Whereas this part of the crossover curve can easily be 
reached by simulating large system sizes with very small interaction ranges, higher-order range dependences prevent 
the direct use of these data for the construction of crossover curves. It was recognized in Ref. || that these are the 
same corrections that are responsible for the deviations from the asymptotic range dependence in Figs. I|, ra , and |(], 
so that this effect can be removed by dividing the magnetization density by the factor in brackets in Eq. (16|) and the 
other quantities by the corresponding counterparts of this factor. 



B. Magnetization density 

As follows from Eq. (^), the magnetization density at criticality is proportional to L Vh ~ 3 in the Ising regime. 

The prefactor depends on the interaction range and scales as R Vh . On the other hand, (|m|) is independent of R 
in the mean-field regime and just scales as TV -1 / 4 oc L~ 3 / 4 . If the crossover behavior can indeed be described in 
terms of a single variable G = L/R 4 , a data collapse should be obtained for (|m|)L 3 / 4 . In the mean-field regime, this 
quantity is independent of G and in the Ising regime it scales as G 2 "*" 9 / 4 . The resulting curve for this quantity is 
shown in Fig. 0. It is immediately clear that the data lie on a perfectly smooth curve, confirming that the crossover 
is indeed ruled by the generalized Ginzburg parameter G. The correction parameter C[m] = 1 + A\R~ 2 refers to the 
higher-order range dependences which have been divided out, in order to make the data for small R collapse on the 
same (Ising) asymptote. For large interaction ranges this correction factor rapidly approaches unity. In the graph 
I have included a line with slope — 9/4 = 0.2315, indicating the dependence on G in the Ising regime. Whereas 
no exact result exists for the finite-size amplitude of this asymptote, it is possible to calculate its counterpart in the 
mean-field regime, where it is found that jg] 

(| m |) i 3/4 = 12 l/4£(|) +0( _i_ ) . (19) 

Thus, (|m|)L 3 / 4 should approach 0.909 891 ... in the limit G — > 0. One indeed observes that the leftmost data points 
in the graph lie already very close to this limit. Together with the collapse of all numerical data onto a single curve, 
this also indicates that the simulations for systems with large interaction ranges indeed have been carried out at the 
correct temperatures; i.e., the extrapolation formula Eq. (|lj) has yielded sufficiently accurate estimates for the critical 
temperatures for 18 < R 2 n < 540. For the sake of clarity, it is stressed that for each single value of R, the simulations 
of the finite systems have been carried out at the critical temperature of a system with that particular interaction 
range in the thermodynamic limit. 

As a side remark, I note that a much more sensitive description of the crossover can be formulated in terms 
of so-called "effective exponents". Originally introduced by Kouvel and Fisher |32|, they have found widespread 
use in experimental analyses (see, e.g., Ref. |24|]) and more recently also in the analysis of numerical results, cf. 
Refs. |7|,p|, ^0p^ ] . Although these effective exponents are usually defined in terms of the logarithmic derivative with 
respect to the reduced temperature, an effective magnetic exponent can be introduced as 

ff _ 9 dH(\m\)L^) 
Vh =4 + d\n(L/R 4 ) ■ (20) 

In the mean-field regime, y£ ff does not approach the classical value y^ = 1 + d/2, but the corresponding value 
y^ = 3d/ 4. This directly related to the violation of hyperscaling in the mean-field regime and can be explained from 
the dangerous- irrelevant- variable mechanism |33"|-|35"|| . This is clearly illustrated in Fig. where a smooth interpolation 
between the value 9/4 and the Ising value 2.4815 is found. 







C. Susceptibility 



In a very similar way, the crossover function for the magnetic susceptibility x a t criticality can be obtained. Since 
it is proportional to the average square magnetization density, it is independent of R in the mean-field regime. In 
the Ising regime, it scales as L 2yh ~ 3 R 2 ( 9 ~ 4:Vh \ so that the quantity xL~ 3 ^ 2 can be represented as a function of the 
parameter G. Indeed, upon application of the range-dependent correction factor C[x], which has the same form as 
the factor between brackets in Eq. (|l6|), a perfect data collapse is obtained, see Fig. [| The total crossover curve 
spans approximately four decades in G, just as for the magnetization density. The exact mean-field result expected 
here is xL~ 3 ^ 2 — ► v^2r(|)/r(~) = 1.170 829 . . ., which is indeed well reproduced for the data in the regime G — ► 0. 
No nonlinear finite-size effects can be observed, suggesting that these are (on the scale of the graph) negligibly small 
for L > 20. 



D. Fourth-order amplitude ratio 



Rather than reproducing crossover curves for the connected susceptibility or the spin-spin correlation function, which 
are very similar to those presented in the previous two subsections, I prefer to pay some attention to the crossover 
of the amplitude ratio Q. This quantity, which is just a disguised form of the fourth-order cumulant introduced by 
Binder [ ^6| , attains trivial limiting values on either side of the critical temperature, but takes a nontrivial universal 
value at criticality. Its Ising limit Q\ — 0.6233 (4) has already played an important role in Sec. IIIB , where this 
parameter was used to determine the location of the critical point. The critical value in the mean- field limit is 
known exactly, Qmf = 0.456 946 58 . . . [f37]j|]. Indeed, the full crossover from Qmf to Q\ as a function of L/R 4 can 
be observed, as illustrated in Fig. [l(]. No correction term has been applied here, because it may be expected that 
the correction terms for (m 2 ) 2 and (m 4 ) cancel each other to a large extent, cf. Fig. 8 in Ref. Q. The less smooth 
appearance of the crossover curve compared to that for the magnetization density and the susceptibility can mainly be 
attributed to several other effects. Apart from the much larger scale of the graph, it turns out that nonlinear finite-size 
effects arc considerably stronger for Q than for other quantities. Further deviations are caused by imperfections in 
the estimates for T c for large R, which on this scale become visible for the larger system sizes. 



VI. CONCLUSIONS 



In this paper, I have presented a detailed determination of the critical properties of the three-dimensional equivalent- 
neighbor model, which is a generalization of the spin-i Ising model, on a cubic lattice. Monte Carlo simulations have 
been carried out for systems with up to thirteen neighbor shells, corresponding to 250 equivalent neighbors. All systems 
have been shown to belong to the 3D Ising universality class. An analysis of these critical properties has yielded a 
coherent picture of their dependence on the interaction range R. The shift of the critical temperature as a function of 
interaction range, to which various mechanisms appear to contribute, has been determined and compared to theoretical 
predictions. I have shown that the range dependence of the critical finite-size amplitudes of the magnetization density 
and the magnetic susceptibility conform very well to the theoretically expected behavior. Also renormalization-group 
predictions for the variation of the finite-size corrections with interaction range have been confirmed and an estimate 
has been obtained for the effective interaction range at which the leading finite-size corrections should vanish. The 
numerical results support the expectation that the (f> 4 coefficient in the Landau-Ginzburg- Wilson Hamiltonian varies 
monotonically with interaction range (or coordination number) and scales for large ranges as l/R 4 . Further Monte 
Carlo results for systems with very large coordination numbers could be obtained by means of an efficient simulation 
scheme. These results enabled the mapping of the full finite-size crossover curves for several quantities, including the 
magnetic susceptibility and the fourth-order amplitude ratio. All these curves can be described by a single crossover 
parameter L/R 4 and interpolate smoothly between mean-field and Ising- like behavior. Also the finite-size crossover 
function for the effective magnetic exponent yh has been obtained. 

A very interesting and experimentally most relevant extension of the work presented here is the case of thermal 
crossover, for which some first results have appeared in Ref. pC| |. A more extensive analysis of this case will be 
presented elsewhere M. 
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FIG. 1. The leading thermal coefficient in the finite-size expansion for the amplitude ratio Q, as a function of coordination 
number. The dashed curve shows the RG prediction (valid in the large-q limit) of Ref. [pj. In order to appreciate the quality 
of this prediction, a linear q dependence is shown as well. 
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FIG. 2. Range dependence of the leading irrelevant field, cf. the second argument in the right-hand side of Eq. (^). Note 
that the result for R — 1 is not shown, because it has the opposite sign. The dashed line represents the asymptotic expression, 
bi oc ]%~ 4y i/( 4 ~ d ') as discussed in the text. 
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FIG. 3. Numerical results for the inverse critical temperature, normalized by the mean-field critical temperature, as a function 
of the inverse squared interaction range, together with the series-expansion results of Dalton and Domb jl4| and Thouless pj . 
The dashed and the dotted lines indicate the results of the least-squares fits discussed in Section IV A, where the dotted line 
is the phenomenological description in which lattice effects have been ignored. 
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FIG. 4. Range dependence of the critical finite-size amplitude of the magnetization density, together with the predicted 
asymptotic range dependence (dashed line) and a fit of all the data points to the renormalization-group prediction (dotted 
curve) . 
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FIG. 5. Range dependence of the critical finite-size amplitude of the magnetic susceptibility, together with the predicted 
asymptotic range dependence (dashed line) and a fit of the data points to the renormalization-group prediction (dotted curve). 
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FIG. 6. Range dependence of the critical finite-size amplitude ro(R) of the connected susceptibility, together with the 
predicted asymptotic range dependence and a fit of the data points to the renormalization-group prediction. Also the fre- 
quently-used scaled susceptibility is shown, which clearly exhibits stronger deviations from the asymptotic range dependence. 
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FIG. 7. Finite-size crossover curve for the absolute magnetization density {\m\) multiplied by an appropriate power of the 
system size. For very small interaction ranges (rightmost data points), higher-order range dependences have been divided out, 
as indicated by the correction factor C[m] (for a more extensive discussion of this topic the reader is referred to the text). The 
crossover curve spans at least four decades in the parameter L/R 4 and systems with a coordination number up to q — 52514 
had to be employed to fully reach the mean-field limit. The perfect collapse of all interaction ranges and system sizes confirms 
the validity of the crossover description in terms of a single parameter. The dashed lines denote the exact mean-field limit 
("MF") and the Ising asymptote with slope yu — 9/4. 
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FIG. 8. The crossover behavior of the effective magnetic exponent as a function of the finite-size crossover parameter. 
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FIG. 9. Finite-size crossover curve for the magnetic susceptibility multiplied by an appropriate power of the system size. For 
very small interaction ranges (rightmost data points), higher-order range dependences have been divided out, as indicated by 
the correction factor C[x\- Just as in Fig. 0, systems with a coordination number up to q — 52514 had to be employed to fully 
reach the mean-field limit. The perfect collapse of all interaction ranges and system sizes confirms the validity of the crossover 
description in terms of a single parameter. The dashed lines denote the exact mean-field limit ( "MF" ) and the Ising asymptote 
with slope 2yh — 9/2. 
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TABLE I. The range of interaction R m , the corresponding number of neighbors q, and the effective range of interaction R 
for the thirteen neighbor shells examined in this work. 



Shell 


r>2 


<7 


R 2 


1 


1 


6 


1 


2 


2 


18 


5 
3 


3 
4 
5 
6 
7 
8 


3 
4 
5 
6 
8 
9 


26 
32 
56 
80 
92 
122 


27 
13 
39 
16 
99 
28 
171 
40 
219 
46 
354 
61 


9 


10 


146 


474 
73 


10 


11 


170 


606 
85 


11 


12 


178 


654 
89 


12 


13 


202 


810 
101 


13 


14 


250 


1146 
125 



TABLE II. The amplitude ratio Q and critical coupling K c for the various ranges of interaction studied in this paper. The 
numbers in parentheses denote the errors in the last decimal places. The results for R? m = 1 (3D nearest-neighbor Ising model) 
stem from Ref. Jilj . The fourth column shows the estimates for K c obtained with Q fixed at the value found in the same work 
(the error margins include the uncertainty in Q). For comparison, the estimates for K c given in Ref. |lj are listed as well. 



R 2 


Q 


K c 


K c 


K c p] 


1 


0.6233 (4) 


0.2216546 (10) 




0.22171 


2 


0.6238 (8) 


0.0644223 (5) 


0.0644220 (5) 


0.06450 


3 


0.6233 (8) 


0.0430381 (4) 


0.0430381 (4) 


0.0432 


4 


0.6224 (5) 


0.03432668 (12) 


0.03432685 (15) 




5 


0.6216 (14) 


0.01892909 (7) 


0.01892915 (4) 




6 


0.621 (3) 


0.01307105 (7) 


0.01307111 (3) 




8 


0.617 (4) 


0.01130202 (8) 


0.01130213 (3) 




9 


0.608 (10) 


0.00844691 (12) 


0.00844703 (4) 




10 


0.614 (11) 


0.00702798 (9) 


0.00702798 (4) 




11 


0.61 (2) 


0.00601661 (14) 


0.00601663 (5) 




12 


0.624 (11) 


0.00574107 (7) 


0.00574110 (4) 




13 


0.618 (8) 


0.00504666 (3) 


0.00504666 (2) 




14 


0.600 (14) 


0.00406419 (4) 


0.00406422 (2) 





TABLE III. The leading correction amplitudes appearing in the Wegner expansion for the magnetization (T < T c , a m ), the 
magnetic susceptibility (T > T c , a x ) and the squared correlation length (T > T c , a??), for three different lattice structures. 
The results for a m were taken from Ref. jS) and the results for a x and a ? 2 from Re f. Egjl . The (slight) nonmonotonicity as 
a function of coordination number in the latter two quantities, already noted in Ref. |22)7is probably not significant and also 
appears to depend on the adopted choice for the susceptibility exponent 7 (the present results correspond to 7 = 1.237). The 
results for a m correspond to the somewhat too high value B — 0.3305, which can probably account for the difference with the 
result a m ~ —0.203 (for the sc lattice) obtained in Ref. |B8j. 





sc (<j = 6) 


bcc (q = 8) 


fee (q = 12) 




-0.256 


-0.240 


-0.234 


a x 


-0.108 


-0.119 


-0.114 


ftc2 


-0.363 


-0.217 


-0.222 
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TABLE IV. The magnetic exponent y^ and the critical amplitude do(R) of the absolute magnetization density as a function 
of interaction range. The estimates for yu in the third column have been obtained with K c fixed at their best values given in 
Table D, whereas the critical amplitudes have been obtained with yu fixed at its 3D Ising value. 



p2 
iX m 


Vh 


Vh 


d {R) 


2 


2.479 (2) 


2.479 (1) 


0.9674 (5) 


3 


2.479 (2) 


2.481 (1) 


0.8933 (6) 


4 


2.475 (5) 


2.479 (1) 


0.8424 (4) 


5 


2.477 (4) 


2.480 (1) 


0.7269 (5) 


6 


2.476 (6) 


2.483 (2) 


0.6716 (7) 


8 


2.472 (7) 


2.484 (3) 


0.6415 (9) 


9 


2.46 (2) 


2.480 (3) 


0.5895 (10) 


10 


2.47 (2) 


2.478 (3) 


0.5622 (10) 


11 


2.47 (2) 


2.471 (5) 


0.5395 (14) 


12 


2.53 (4) 


2.485 (6) 


0.5335 (20) 


13 


2.47 (2) 


2.480 (5) 


0.5128 (16) 


14 


2.463 (15) 


2.475 (4) 


0.4845 (17) 



TABLE V. The magnetic exponent yu and the critical amplitude po{R) of the magnetic susceptibility as a function of 
interaction range. The estimates for y^ in the third column have been obtained with K c fixed at their best values given in 
Table [ll|, whereas the critical amplitudes have been obtained with yn fixed at its 3D Ising value. The data point for R^ — 1 is 
taken from Ref. [Ell. 



R 2 


Vh 


Vh 


Po(R) 


1 






1.5580 (15) 


2 


2.479 (1) 


2.479 (1) 


1.1620 (7) 


3 


2.481 (6) 


2.484 (3) 


0.9865 (32) 


4 


2.478 (6) 


2.484 (2) 


0.8752 (12) 


5 


2.481 (8) 


2.481 (3) 


0.6518 (18) 


6 


2.478 (13) 


2.478 (12) 


0.5534 (35) 


8 


2.484 (14) 


2.480 (2) 


0.5105 (16) 


9 


2.46 (3) 


2.476 (9) 


0.4343 (12) 


10 


2.46 (2) 


2.474 (4) 


0.3951 (15) 


11 


2.46 (2) 


2.47 (1) 


0.3653 (16) 


12 


2.48 (2) 


2.481 (6) 


0.3564 (24) 


13 


2.46 (2) 


2.484 (5) 


0.3297 (16) 


14 


2.45 (4) 


2.477 (6) 


0.2943 (23) 
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